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The ferrimagnetic spin-1/2 chain composed of alternating Ising and Heisenberg spins in an arbitrarily ori- 
ented magnetic field is exactly solved using the spin-rotation transformation and the transfer-matrix method. 
It is shown that the low-temperature magnetization process depends basically on a spatial orientation of the 
magnetic field. A sharp stepwise magnetization curve with a marked intermediate plateau, which emerges 
for the magnetic field applied along the easy-axis direction of the Ising spins, becomes smoother and the 
intermediate plateau shrinks if the external field is tilted from the easy-axis direction. The magnetization curve 
of a polycrystalline system is also calculated by performing powder averaging of the derived magnetization 
formula. The proposed spin-chain model brings an insight into high-field magnetization data of 3d-4f bimetallic 
polymeric compound Dy(N03)(DMSO)2Cu(opba)(DMSO)2, which affords an interesting experimental realiza- 
tion of the ferrimagnetic chain composed of two different but regularly alternating spin-1/2 magnetic ions Dy 3 + 
and Cu 2 + that are reasonably approximated by the notion of Ising and Heisenberg spins, respectively. 
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1. Introduction 

Exactly solved quantum spin chains belong to the most attracting issues to deal with in the 
condensed matter theory, because they are capable of providing a deeper understanding into many 
unconvential quantum cooperative phenomena Recently, the particular research interest has 
been devoted to sophisticated Ising-Heisenberg chains, which are aimed at describing hybrid spin 
systems composed of 'classical' Ising and quantum Heisenberg spins [^-flol| . Among other matters, 
the Ising-Heisenberg chains have become helpful in order to provide evidence of several novel and 
unexpected quantum states 043 > fractional magnetization plateaus in the low-temperature mag- 
netization process [^-Isl, enhanced magnetocaloric effect during the adiabatic demagnetization [5( , 
thermal entanglement [7|, etc. It is therefore quite challenging to search for suitable experimental 
realizations of the Ising-Heisenberg chains testifying the aforementioned theoretical findings, but 
only a few experimental systems satisfy a very specific requirement of a regular alternation of 
the Ising and Heisenberg spins. Up to now, the magnetic behaviour of three different polymeric 
chains Cu(3-Clpy) 2 (N 3 ) 2 @|, [(CuL) 2 Dy][Mo(CN) 8 ] Q and [Fe(H 2 0)(L)][Nb(CN) 8 ][Fe(L)] [lQj was 
successfully interpreted within the framework of the Ising-Heisenberg chains. 

The main goal of the present work is to examine the magnetization process in the spin-^ chain 
of alternating Ising and Heisenberg spins, which brings an insight into ferrimagnetism of 3d-4f 
bimetallic coordination compound Dy(N0 3 )(DMSO) 2 Cu(opba)(DMSO) 2 The organization of 
this paper is as follows. Exact results for the total and sublattice magnetizations of the investi- 
gated spin-chain model are derived in Section [3] The most interesting theoretical results are then 
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presented in Section [3j where they are also confronted with the relevant experimental data. The 
paper ends up with several concluding remarks given in Section 0J 

2. Spin alternating chain 

Let us consider the spin-^ chain composed of alternating Ising and Heisenberg spins in an 
external magnetic field of arbitrary spatial direction. It is quite plausible to suppose that magnetic 
properties of this 'classical-quantum' spin-chain model will be highly anisotropic due to a presence 
of the Ising spins. In this respect, it is of particular interest to examine how the magnetization 
process depends on a spatial orientation of the external magnetic field, which can be unambiguously 
given by the deviation angle referred with respect to a unique easy axis of the Ising spins. The 
investigated spin-chain model can be defined through the following Hamiltonian 

N N N N 

H = -JY / S* K + <+i ) - 9tmB cos 9 J2 < - 92 MbS sin £ Sf - g 7 2 ^B cos 0^5?. (1) 

i— 1 i— 1 i—l i— 1 

Here, of and Sf (a = x, z) denote standard spatial components of the spin-^ operator, whereas 
the former (latter) operators apparently refer to the Ising (Heisenberg) spins. The first summation 
takes into account the Ising-type exchange interaction J between the nearest-neighbour Heisenberg 
and Ising spins, the second term determines the Zeeman's energy of the Ising spins in the external 
magnetic field with the projection B cos 9 towards their easy (z) axis, while the last two Zeeman's 
terms determine the overall magnetostatic energy of the Heisenberg spins affected both by the 
transverse (Bs'mO) as well as longitudinal (BcosO) component of the external magnetic field. 
The quantities gf and g% (a = x, z) label spatial components of Lande ^-factors of the Ising and 
Heisenberg spins, respectively, is Bohr magneton and B stands for the external magnetic field. 
Notice that the Hamiltonian (TTJ) is built on assumption that the transverse component of Lande 
g- factor of the Ising spins is negligible (gf « 0), i.e. the situation, which is realistic only for the 
highly anisotropic (the so-called Ising-type) magnetic ions (l2l . Il3| . 

Taking advantage of a 'classical' nature of the Ising spins, which represent barrier for local quan- 
tum fluctuations induced by the transverse component of the external field acting on the Heisenberg 
spins, one may rewrite the total Hamiltonian (fT]) as a sum of commuting site Hamiltonians 

N 

U = YKi, (2) 

whereas each site Hamiltonians Hi involves all the interaction terms of the ith Heisenberg spin 
and the Zeeman's energy of its two nearest-neighbour Ising spins 

n t = -JSl(ol + a? +1 ) -H%S*-HlS*-?± « + . (3) 

For simplicity, we have introduced here three 'effective fields' Hf — gl^B cos 9, fff = g%[i^,B sin 9, 
and -fff = gf^B^cosf? in order to write the Hamiltonian and all subsequent expressions in a 
more abbreviated form. Owing to a validity of the commutation relation [Hi, Hj] = between dif- 
ferent site Hamiltonians, the partition function of the considered spin-chain model can be partially 
factorized into the following product 

N N 

z = e n Tr ^ ex p (-w = e n t ( 4 ) 

{a,}i=l {o-«}i=l 

where (3 = l/(k-B,T). fee is Boltzmann's constant, T is the absolute temperature, the symbol Trs ; 
denotes a trace over two spin states of the ith Heisenberg spin and the summation } runs over 
all available configurations of the Ising spins. To proceed further with a calculation, the partial 
trace over spin degrees of freedom of the Heisenberg spins must be performed before summing over 
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spin states of the Ising spins. It is therefore quite convenient to diagonalize the site Hamiltonian 
([3]) by making use of the spin-rotation transformation 



Sf = Sf cos ^ + Sf sin (/),, Sf = -Sf sin < 
which brings the site Hamiltonian ([3]) into the diagonal form 



Sf cos < 



-Hi = -5fV [J(a| + af +1 ) + H^ 2 + (iff )a - ( CTf + 
provided that the spin rotation is performed by the specific angle 

m 

— arctan 



(5) 



(6) 



(7) 



The effective Boltzmann's weight, which enters the factorized form of the partition function 
([!]), can be now simply evaluated by employing the trace invariance and the diagonalized form of 
the site Hamiltonian ^ 



T(of,of +1 ) = Tr Si exp(-/^) = Trs!exp (_^/) 



(8) 



2 exp 



cosh 1 1 ^[JK +< +1 ) + H|] 2 + (i/f) 



The effective Boltzmann's factor ([8]) apparently depends, after tracing out the spin degrees of 
freedom of the ith Heisenberg spin, only upon its two nearest-neighbour Ising spins Oi and (7j.fi • 
Thus, the expression ([5]) can alternatively be viewed as the effective two-by-two transfer matrix 



2' 1 2 
1 i 1 



2' 1 2 

with three different matrix elements defined as 



l __L 

'I' 2 

I _I 

2 ' 2 



T+ To 
T T_ 



(9) 



1 , 1 



T ± ee T ±-, ±- = 2 exp ±^fff cosh 



/3 



^/(tff) 2 + (tff) 2 



(10) 



Substituting the matrix ^ into the relation ^ one may consequently employ the standard 
transfer-matrix approach in order to obtain the exact result for the partition function 



N 



z =En T «'< + i)= TrTAr=A ^+ A 



N 
2 I 



which is written in terms of two respective eigenvalues of the transfer matrix (JSJ) 

Al.2 = 



TV 



-4T 2 



(11) 



(12) 



Now, let us proceed to a calculation of the most important quantities, which are relevant for our 
subsequent analysis of the magnetization process. The Gibbs free energy can easily be calculated 
from the exact expression pip for the partition function. In the thermodynamic limit N —¥ oo, 
one obtains the following precise analytical result for the free energy per elementary unit 



f = —k^T lim — InZ = A;BTln2 — /cbTIii 



T+ + T_ + J(T+ - T_) 2 + 4T 2 



(13) 



Subsequently, one may readily calculate the sublattice magnetizations of the Ising and Heisenberg 
spins by differentiating the free energy (I13[) with respect to the appropriate effective fields. The 
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sublattice magnetization of the Ising spins in two mutually orthogonal directions of the external 
magnetic field oriented either perpendicular or parallel with respect to the easy axis read 



mf = 5i>B«) =0, ml = glm(<? z i) 



T, 



2 v/(r + -T_ 



4T 2 



(14) 



Similarly, the sublattice magnetization of the Heisenberg spins in two aforementioned orthogonal 
directions of the external magnetic field can be calculated from the following unique formula 



52 >b (T+ - T_)(t7? - U°) + 4T (7 Q + (U% + U«)^{T+-T_) 



4T 2 



4T 2 



(15) 



2 (T+ - T_) 2 + 4T 2 + (T+ + T_) ^(T+ - T_) 2 

which is valid for both spatial directions a — x and z if the coefficients U± and Z/q* are defined as 

mi 



TT X 



HI 



2 exp ± 



sinh 



(J±fff) 2 



z 2 sinh 



(my 



and 



Ul = ±- 



t/rT 



HI 



(Hf)2 



(F|) 2 
1 2 sinh 



2 exp ± 



(i/f) 2 



sinh 



(16) 



(J±H|) 2 + (i/f) 2 



my 



(17) 



It should be pointed out that the exact analytical formula for the free energy (fT5)) allows a 
straightforward derivation of the magnetization formula for the most general case of arbitrarily 
oriented external magnetic field as well. The final magnetization formula in the external field of 
arbitrary spatial direction can be conveniently expressed in terms of the formerly derived sublattice 
magnetizations (fl"4"|) and (|15[) of the Ising and Heisenberg spins 

m(9) = (m\ + mjj) cos6> + vn\ sin0. (18) 

It is noteworthy that the magnetization formula (|18p might be of fundamental importance for an 
analysis of the angular dependence of the magnetization process in a single-crystal sample related 
to the considered spin-chain model. On the other hand, the magnetization of a powder sample 
pertinent to the investigated spin-chain model can also be formally obtained by integrating the 
magnetization formula (|18p over one hemisphere yielding 



ID V 



m(9) sinfldfl, 



(19) 



but the involved integral precludes derivation of the closed-form analytical expression because of 
too complicated functions involved in the sublattice magnetizations (|14p and (|15p of the Ising and 
Heisenberg spins, respectively. Of course, the integral entering the relation (fill)) for the magnetiza- 
tion of powder samples can be evaluated numerically and hence, it may be of practical importance 
for an investigation of the magnetization process in the related polycrystalline systems. 



3. Results and discussion 



Let us proceed to a discussion of the most interesting findings acquired for the ferrimagnetic 
spin-i chain of alternating Ising and Heisenberg spins coupled through the antiferromagnetic 
nearest-neighbour interaction J = —\J\ < 0. First, we will present a comprehensive survey of 
theoretical results with the aim to shed light on typical magnetization features of the proposed 
spin-chain model and then, high-field magnetization data of one specific polymeric coordination 
compound will be clarified within the framework of the model under investigation. 
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3.1. Survey of theoretical results 



Let us begin with the analysis of the ground state. The diagonal form of the site Hamiltonian 
(JHJ) allows one to get the lowest-energy eigenstate of the investigated spin alternating chain 



N i 1 

ifri) = n k = - 



Sf = - 



Sf = 



(20) 



which indicates the quantum ferrimagnetic ordering with a perfect alignment of the Ising spins 
towards their easy axis and respectively, the quantum superposition of two spin states of each 
Heisenberg spin that basically depends on a mutual interplay between the exchange constant, 
Lande g-factor, size and spatial orientation of the external field via the occurrence probabilities 



it 



|J| — <7f/jB-Bcos# 



\J {\A ~~ gfMB-Bcosfl) 2 T (gf^B-Bsin#)^ 



(21) 



The perfect alignment of the Ising spins is also confirmed by the following zero-temperature values 
of the relevant sublattice magnetization in two conspicuous orthogonal directions (parallel and 
perpendicular) with respect to the easy axis 



« 1 . 
m i = 20iMb, 



0. 



(22) 



Contrary to this, both spatial components of the sublattice magnetization of the Heisenberg spins 
are subject to the quantum reduction of magnetization on behalf of local quantum fluctuations 
arising from the transverse field (i.e. perpendicular projection of the external magnetic field with 
respect to the easy axis of the Ising spins) 



rn 9 = — - 



92 m 



\J\ - sf/XB-Bcos0 



\J\ — .gf/jB-B cos#) 2 + (g% MB-Bsin#) 2 



gf/jB-Bsinfl 



(| J| — gf/iB-Bcos^) 2 T (gf ^b-B sm9) 2 



(23) 



It is quite obvious from (|2"2"|) and (|2"3"|) [or alternatively from (|2"0|) and (|2"Tj) ] that the quantum 
ferrimagnetic order changes to the classical ferrimagnetic order whenever the transverse projection 
of the external magnetic field vanishes. However, the total magnetization of the ferrimagnetic chain 
of alternating Ising and Heisenberg spins is subject to the quantum reduction of magnetization for 
any other spatial orientation of the external magnetic field 



m 



/jb cos 9 



32Mb sin( 



91 - 92 



\J\ — gf /ib-B cos 6 



(\J\-S5HbB cos9) 
g'zHBB sin 6 



(gf Mb-B sin( 



| J| — <7|/ibB cos0) 2 + (gf/JB-Bsin^) 2 



(24) 



which appears due to the quantum entanglement of two spin states of each Heisenberg spin aris- 
ing from the transverse field. In accordance with this statement, the quantum reduction of total 
magnetization becomes the greater, the higher the transverse component of the external field is. 

Now, let us turn our attention to a discussion of the low-temperature magnetization process 
serving in evidence of the aforedescribed ground-state features. For simplicity, we will further 
assume that the Heisenberg spins have the completely isotropic Lande factor gf = g| = 9% = 2 in 
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Figure 1. Longitudinal and transverse projections of the sublattice magnetization of the Ising 
(ml, mf = 0) and Heisenberg (m|, rn%) spins as a function of the magnetic field at low enough 
temperature ksT/\J\ = 0.01, one particular choice of p-factors gf = 6, gf = 0, <?f = sf = 2, and 
several spatial orientations of the external magnetic field. 



contrast to the highly anisotropic Lande factor of Ising spins gf 3> 2 and gf = g\ = 0. To provide 
in-depth understanding of the magnetization process, figure Q] illustrates typical field dependences 
for two orthogonal projections of the sublattice magnetization of the Ising and Heisenberg spins 
under different spatial orientation of the applied magnetic field. If the external field is applied 
along the easy axis of the Ising spins, one observes a steep field-induced increase in the longitudinal 
sublattice magnetization of the Heisenberg spins m| at a critical field due to an abrupt reversal of all 
Heisenberg spins towards the external-field direction (at low fields, the Heisenberg spins are aligned 
antiparallel with respect to the external field, because they possess lower ^-factor compared to that 
of the Ising spins). This classical mechanism for a formation of the intermediate magnetization 
plateau will consequently lead to a sharp stepwise profile in the total magnetization versus the 
external-field dependence (see figure [5] and the subsequent discussion). If the external field is tilted 
from the easy-axis direction, however, the relevant behaviour of the Heisenberg spins becomes 
much more intricate owing to a mutual competition between two different spatial directions given 
by the easy axis of the Ising spins and the spatial orientation of the applied magnetic field. As a 
matter of fact, the longitudinal component of the sublattice magnetization of the Heisenberg spins 
m| apparently exhibits a smoother variation if the external field is gradually tilted from the easy 
axis and the more gentle field-induced reversal of the Heisenberg spins (given by the change of sign 
of mf ) simultaneously appears at the higher critical field 



1 



Mb -Be 

|J| gl cos 8' 



(25) 



It should be mentioned, moreover, that the transverse projection of the sublattice magnetization 
of the Heisenberg spins mf becomes non-zero for this more general case and mf exhibits a striking 
non-monotonous dependence with a more or less sharper global maximum located at the critical 
field (|2~5|) . The particular case of the external magnetic field oriented perpendicular with respect 
to the easy axis deserves a special mention. Under this specific constraint, the Ising spins become 
completely uncorrelated (disordered) and the Heisenberg spins gradually tend to align towards the 
external-field direction upon strengthening the transverse field. 
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Figure 2. Low-temperature dependence of the total magnetization as a function of the magnetic 
field for several spatial orientations of the external field and two different sets of g-factors: (a) 
gf = 6, gf = 0, gf = 92 = 2; (b) gf = 20, gf = 0, gf = gi = 2. Broken lines show magnetization 
curves for the powder samples, which were obtained by the numerical solution of (119|) , 

In figure [5J the total magnetization is plotted against the magnetic field for various spatial 
orientations of the external field and two different values of the Lande g-factor of the Ising spins. 
It is worthy of notice that the profile of displayed magnetization curves can readily be understood 
from the relevant low-temperature dependences of the sublattice magnetizations (cf. figure [2] with 
figure [T]). In agreement with our expectations, the classical mechanism for a formation of the 
intermediate magnetization plateau becomes quite evident under the particular orientation of the 
external field directed along the easy axis of the Ising spins. Under this specific condition, one 
actually detects a sharp stepwise magnetization curve with a marked intermediate plateau at the 
fractional value (gf — g|)^B/2 per elementary unit, which implies an existence of the classical 
ferrimagnetic order due to a unequal magnetic moment of the nearest- neighbour spin-^ atoms 
differing in their <?-factors. It can be clearly seen from figure[5]that the intermediate magnetization 
plateau gradually shrinks and the relevant dependence of the total magnetization becomes smoother 
as the external field deviates from the easy-axis direction of the Ising spins. The observed breakdown 
of intermediate magnetization plateau, the gradual smoothing of the magnetization curve as well 
as the overall quantum reduction of the total magnetization can all be attributed to local quantum 
fluctuations, which arise from the transverse component of the external magnetic field acting on the 
Heisenberg spins. For comparison, the low-temperature magnetization curve of a polycrystallinc 
system is also depicted in figure [5] by a broken line. Interestingly, the powder averaging through the 
formula (|19p yields the magnetization curve of a polycrystalline system, which quite closely follows 
the magnetization curve of a single-crystal system for the particular orientation of the external 
magnetic field deviating by the angle 9 — 60° from the easy-axis direction. 

3.2. High-field magnetization of DyCu 

Our theoretical findings for the magnetization process will be now confronted with high-field 
magnetization data of 3d-4f bimetallic coordination polymer Dy(N03)(DMSO)2Cu(opba)(DMSO)2 
(DMSO=dimethylsulfoxide, opba=orthophenylenebisoxamato) to be further abbreviated as DyCu. 
The polycrystalline sample of DyCu was prepared according to the procedure reported on previ- 
ously by Calvez and co-workers Even though our effort aimed at a preparation of a single- 
crystal sample suitable for a crystal-structure characterization and magnetization measurements 
was not successful, the elemental analysis has provided a strong support that the prepared polycrys- 
talline sample ought to consists of bimetallic polymeric chains representing a structural analogue 
of Ln(N0 3 )(DMSO) 2 Cu(opba)(DMSO) 2 (Ln=Gd-Er), which is visualized in figure H by making 
use of the crystallographic data reported by Calvez et al. [IlJ . Bearing this in mind, it could be ex- 
pected that the bis-chelating bridging ligand opba might mediate a relatively strong superexchange 
coupling between the nearest-neighbour Dy 3+ and Cu 2+ magnetic ions, which should consequently 
form the bimetallic polymeric chain running along the crystallographic c-axis. 
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Figure 3. The visualization of LnCu polymeric chain in bimetallic coordination compounds 
Ln(N0 3 )(DMSO) 2 Cu(opba)(DMSO)2 (Ln=Gd-Er) by adopting the crystallographic data from 
Ref. [lH deposited at The Cambridge Crystallographic Data Centre. Coloring scheme for the 
atom labelling: Ln (magenta), Cu (green), O (red), N (blue), C (grey), S (yellow). 
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Figure 4. (a) High-field magnetization data recorded for the powder sample of DyCu at two 
sufficiently low temperatures 1.3 K and 4.2 K; (b) High-field magnetization curve of DyCu at 
the lowest measured temperature 1.3 K before and after subtracting the paramagnetic Van Vleck 
contribution estimated from the quasi-linear dependence in the high-field range 36-52 T. 



First, let us make a few comments on constituent magnetic ions of DyCu. It is quite well 
established that the magnetic behaviour of octahedxally coordinated Cu 2+ ions can be quite faith- 
fully represented by the notion of Heisenberg spins with a more or less isotropic Lande factor 
g% g% « <j>! > 2 [13], while the eight-coordinated Dy 3+ ions often obey more subtle requirements 
of the Ising spins necessitating a highly anisotropic ^-factor g\ 3> 2 and gf « gf sa @, 12- H]. 



In fact, Dy 3+ ion represents Kramers ion with the ground-state multiplet 6 H 15 / 2 , which us ually 
undergoes a rather strong crystal-field splitting into eight well-separated Kramers doublets [l2lll3j. 
In this respect, the magnetic behaviour of Dy 3+ ion can be often interpreted at low enough temper- 
atures in terms of the Ising spin with the effective Lande factor gf s» 20 and gf « g\ f=s (JJ ll2HT4 1 . 
With all this in mind, the coordination polymer DyCu could thus provide a suitable experimental 
realization of the spin-^ chain of the alternating Ising and Heisenberg spins. However, one should 
also keep in mind that this simplification is reasonable only at sufficiently low temperatures com- 
pared with the energy gap between the lowest-energy and excited Kramers doublets to avoid any 
danger of over-interpretation inherent in this approximation. 

Figure @Ia) shows the high- field magnetization curve of the powder sample of DyCu measured 
at two sufficiently low temperatures 1.3 K and 4.2 K (see reference (TBJ for more details on the 
set-up of magnetization measurements). The displayed magnetization curves exhibit a remarkable 
crossing around 12 T, whereas the magnetization data recorded at a higher temperature 4.2 K 
become greater than the ones recorded at a lower temperature 1.3 K above this crossing field. 
Apart from this fact, one may clearly recognize three characteristic regions in the relevant magne- 
tization process as illustrated in figure B|b) by thin solid lines serving as guides for eyes only. The 
magnetization initially shows a very rapid increase with the magnetic field in the low-field range 
0-3 T, then it exhibits a rather steep increase in the range of moderate fields 5-32 T, which is 
consecutively followed with a more steady quasi-linear dependence in the high-field range 34-52 T. 
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(a) S[T] (b) B[T] 



Figure 5. High-field magnetization data of the powder sample of DyCu at two different tem- 
peratures: (a) T = 1.3 K; (b) T = 4.2 K. Blue solid lines show the best theoretical fit obtained 
using the formulas (|18[) and (|19|) for the spin-^ chain of alternating Ising and Heisenberg spins. 




Figure 6. High-field magnetization data of the powder sample of DyCu at two different tem- 
peratures: (a) T — 1.3 K; (b) T = 4.2 K. Blue solid lines show the best theoretical fit obtained 
using the formulas (|26l) and (j 19ft for the spin-^ chain of alternating Ising and Heisenberg spins 
with the mean-field correction for the transverse magnetization of the Ising spins. 



The steady quasi-linear dependence observed in the high-field range implies a significant contribu- 
tion of the temperature-independent (Van Vleck) paramagnetism [3] , which was evaluated to be 
wtip = 0.029/xb.T -1 per Dy 3+ ion from the linear fit of the high-field region 36-52 T and subse- 
quently subtracted from the measured magnetization data. Even though this procedure might give 
only a rather rough estimate of mxip, the obtained value is in a rather good accord with typical 
values of totip reported on previously for other compounds involving Dy 3+ ion Q. 

High-field magnetization data of the powder sample of DyCu after subtracting the temperature- 
independent paramagnetism are presented in figure [5] together with the best theoretical fit obtained 
by using the magnetization formulas (|18[) and (|19j) derived for the ferrimagnetic spin-^ chain of 
alternating Ising and Heisenberg spins. As one can see, one generally obtains a quite plausible 
concordance between the recorded experimental data and the relevant theoretical predictions ex- 
cept a small discrepancy observable in the low-field region, where a more abrupt change in the 
magnetization is theoretically predicted than experimentally observed. It should be nevertheless 
mentioned that the determined value of Lande factor of Dy 3+ ion g\ = 22.6 is somewhat greater 
than theoretically expected. This might indicate a small but non-negligible effect of the transverse 
component of (/-factor , which was completely ignored for Dy 3+ ions within our exact treatment. If 
the transverse component gf of Lande factor is taken into account for Dy 3+ ions at the mean-field 
level, the magnetization formula (jT5J) for a single-crystal sample can easily be corrected for the 
contribution of the transverse magnetization of the Ising spins (Dy 3+ ions) yielding 

m(9) = (mf + ml) cos + (mf + mf) sin 6 (26) 
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with 



y/(2Jm% + H{ ) 2 + (iff ) 2 2 



tanh 



(27) 



In doing so, the corrected magnetization formula (|26p can be substituted into the formula ([TT))) 
derived for the polycrystalline system in order to analyze the magnetization curve of the powder 
sample of DyCu. In this way, one actually resolves the problem with too high value of the Lande 
g-factor of Dy 3+ ions as evidenced by the best fitting data set indicated in figure [51 It is also quite 
evident from figure [5] that the acquired theoretical prediction follows more accurately the relevant 
experimental data in the low-field region and beside this, the values of ^-factors gained for Dy 3+ 
and Cu 2+ ions from this latter fitting procedure are simultaneously much more reliable. 



4. Concluding remarks 

The present article is devoted to an exact study of the spin-i chain of alternating Ising and 
Heisenberg spins in the magnetic field of arbitrary spatial direction. The low-temperature mag- 
netization curve of the investigated spin-chain model was detailed examined in dependence on a 
spatial orientation of the applied magnetic field. It has been demonstrated that the magnetization 
curve becomes smoother, the intermediate plateau shrinks, and the total magnetization is reduced 
by quantum fluctuations as the external field deviates from the easy axis of the Ising spins. Accord- 
ingly, the powder averaging in the related polycrystalline system yields a smooth low-temperature 
magnetization curve, which quite closely follows the magnetization curve of a single-crystal system 
for a spatial orientation of the external field deviating by 9 = 60° from the easy-axis direction. 

The proposed spin-chain model has been employed for interpreting high-field magnetization 
data of polymeric coordination compound DyCu, which affords an interesting experimental real- 
ization of the ferrimagnetic spin-^ chain of regularly alternating Dy 3+ and Cu 2+ magnetic ions 
reasonably approximated by the notion of Ising and Heisenberg spins, respectively. From this 
perspective, experimental studies performed on single-crystal and field-aligned samples of DyCu 
represent a particularly challenging problem for future investigations. 
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